Thermally activated magnetization reversal in monoatomic magnetic chains on 
surfaces studied by classical atomistic spin-dynamics simulations 
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We analyze the spontaneous magnetization reversal of supported monoatomic chains of finite 
length due to thermal fluctuations via atomistic spin-dynamics simulations. Our approach is based 
on the integration of the Landau-Lifshitz equation of motion of a classical spin Hamiltonian at the 
presence of stochastic forces. The associated magnetization lifetime is found to obey an Arrhenius 
law with an activation barrier equal to the domain wall energy in the chain. For chains longer than 
one domain-wall width, the reversal is initiated by nucleation of a reversed magnetization domain 
primarily at the chain edge followed by a subsequent propagation of the domain wall to the other 
edge in a random-walk fashion. This results in a linear dependence of the lifetime on the chain 
length, if the magnetization correlation length is not exceeded. We studied chains of uniaxial and 
tri-axial anisotropy and found that a tri-axial anisotropy leads to a reduction of the magnetization 
lifetime due to a higher reversal attempt rate, even though the activation barrier is not changed. 
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I. INTRODUCTION 

Advances in experimental techniques in recent years 
have made possible the controlled growth and character- 
ization of magnetic chains on non-magnetic crystal sur- 
faces. Such chains can be made as thin as monoatomic, 
i.e., their cross-section consists of a single atom, while 
they are usually a few nanometers or tens of nanome- 
ters long; they grow, e.g., at terrace step edges, in the 
"trenches" of (110) surfaces or as inclusions in surface 
alloys 0-01, so that their structure is rather stable, and 
they can be studied by a number of spin-sensitive tech- 
niques including spin-polarized scanning tunneling mi- 
croscopy or x-ray magnetic circular dichroism. Magnetic 
chains, as all magnetic nanostructures, bear technolog- 
ical relevance due to the prospect of miniaturization of 
magnetic bits for information storage. Particularly ap- 
pealing, however, in chains is that their geometry is de- 
fined by only one parameter, the length, which makes it 
easier to interpret their physical properties @. 

At such small sizes magnetic states at low-temperature 
equilibrium consist of a single domain and the possibil- 
ity of information storage relies solely on the bistability 
of the magnetization, which is a consequence of mag- 
netocrystalline anisotropy. At low enough temperature 
the direction of magnetization is trapped close to a local 
energy minimum for some time before it is reversed by 
thermal flutuations to another minimum, pointing along 
the opposite direction. Above some characteristic block- 
ing temperature Tr, however, there is a transition to a 
situation where the magnetization fluctuations are too 
intense for a stable state to be formed. As the size of 
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the nanomagnets is below the thermodynamic limit, the 
transition is not a genuine phase transition but rather a 
smooth crossover. 

These effects are strongly size dependent. In many 
cases, larger system size leads to higher blocking temper- 
ature, reaching the Curie temperature in the thermody- 
namic limit. However, in two-dimensional [Io| and even 
more in one-dimensional systems, upon increasing the 
system size L, the entropy leads the magnetic state to a 
multi-domain structure at even low temperatures. The 
characteristic length, at which this effect sets in, is the 
magnetization correlation length £(T); thus, the effect of 
magnetic bistability in chains is physically limited by the 
requirement L < £(1b). 

The temperature-driven magnetization reversal dy- 
namics of nano-particles has been tradidionally stud- 
ied within the Stoner-Wohlfarth model [IJ Ell > where 
the particle magnetization is approximated by a single 
"macro-spin" of rigid magnitude. However, it is known 
that excited states departing from the macro-spin ap- 
proximation can play a significant role even at sizes of 
a few nanometers [13[ , while another mechanism sets in, 
namely reversal through domain nucleation and domain 
wall propagation [l4j ; the latter mechanism has been seen 
also in elongated particles PH , and should be even more 
relevant in chains flrjj . 

Spin dynamics simulations of magnetic nanostructures 
are based on a parametrization of the spin-dependent 
part of the total energy and usually follow one of two 
paths: either that of a quantum Hcisenberg model or the 
one of a classical Heisenberg model. It is well-known that 
the quantum model exhausts computational resources at 
relatively small system sizes, as the Hilbert space in- 
creases exponentially with the number of atoms, there- 
fore for system sizes of 50 or 100 atoms only the clas- 
sical model is practically available. However, there are 
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also physical situations where the quantum Hcisenberg 
model is not applicable, in particular magnetic systems 
on metallic surfaces. The quantum model is based on 
the assumption of half-integer or integer spins that in- 
teract with each other. This is the case with spins on 
an insulating surface, such as Mn chains on CuN [l7j |. 
On the other hand, magnetic atoms on metallic surfaces 
are characterized by non-integer or non-half-intcger spin 
values, which are actually average values of the spin op- 
erator on superposition states of the magnetic atoms 
with the substrate continuous spectrum. In this pic- 
ture, where substrate electrons hop on and off the mag- 
netic atoms, the spin states are blurred out compared to 
the sharp eigenstates of the quantum Heisenberg model. 
In addition, damping of the magnetic excitations into 
electron-hole pairs emerges, with concequences in g- shift 
and frequency-dependent magnetization lifetime [18l - l22| . 
This, together with temperature effects, causes decoher- 
ence of the magnetic states and thus the quantum Heisen- 
berg model is not applicable any more; possibly a clas- 
sical model, founded on the principles of adiabatic spin 
dynamics [2J], is better suited in this case. Actually the 
classical atomistic spin dynamics is based on the Landau- 
Lifshitz equations of motion, with the point-dependent 
magnetization being replaced by atom-dependent mag- 
netic moments, and is a widely used tool for the study of 
atom-dependent nanomagnetism [24|, [25[ . 

Having this in mind, we use a classical spin-dynamics 
simulation method, implemented in our juSpinx code, 
to study the time evolution of the magnetization in 
monoatomic chains. Assuming a singlc-ion anisotropy 
tensor Ki, inter- atomic exchange coupling constants Jij 
between atomic spin moments Mj and Mj , and an exter- 
nal magnetic field B cyit , the spin-dependent excitation en- 
ergy of the system is parametrized by a classical Heisen- 
berg Hamiltonian of the form 

H = - - JijMi ■ Mj - Yl ^ h ' K ^ h ~ ^° Xt ' 

ij i i 

(1) 

where it is implied that J a — 0, while the factor ^ com- 
pensates for the double-counting in the summation; i 
and j denote the atoms in the chain. We set henceforth 
\Mi\ = 1, so that Mi are to be understood as unit vec- 
tors along the atomic moment directions, with the mod- 
ulus of the moments absorbed in the parameters Jij , K i 
and £? oxt . Dipole-dipole interactions are neglected here, 
as we are interested in systems of not more than 100 
atoms in size. In order to ensure bistability, the single- 
ion anisotropy should favour a particular axis; the sim- 
plest such case is that Ki is the same for each atom and 
its eigenvalues are (K x ,K y , K z ), with K x = K y = and 
K z = K > 0, i.e. the anisotropy is uniaxial, so that sec- 
ond term of the Hamiltonian takes the form — J^i ^Mf. z , 
with the z-azis appropriately chosen along the low-energy 
direction. In the limit Kj\ | 3> 1, there can be a further 
simplification, as the system behaviour is well-described 
by the Ising model and it follows Glauber dynamics (2|| . 



This condition is met e.g. in single chain magnets [27[ 
consisting usually of transition metal atoms that are sta- 
bilized and linked in a chain form by molecular ligands. 
However, for metallic magnetic chains on surfaces the 
magnetic anisotropy can take large values that are of the 
order of a few meV, while the exchange couplin g is typ- 
iclaly one or two orders of magnitude stronger [28| . In 
this case, an appropriate classical approximation to the 
spin dynamics is the stochastic Landau-Lifshitz equation 
of motion: 

where we have defined the effective magnetic field acting 
on the spin at the atomic site i as 



Bf a = -■ 



dH 



dMi 



(3) 



+ 2 \Ki, x {Mi -x)x + K l , y (M l -y)y + K hZ (M z ■ z) zj , 

and A is a damping parameter leading the system to equi- 
librium at zero temperature. The temperature T enters 
via a stochastic fluctuating force fi(t), as described in 
Appendix A. In Eq. @, the first term on the right-hand 
side is derived directly from the Hamiltonian (p}, while 
the damping and temperature terms are phenomenologi- 
cally added assuming interaction of the spins with a tem- 
perature bath. 

Our focus in this work is on the temperature-induced 
magnetization reversal of finite-length chains. It is im- 
plied that the length of the chain, L, is considerably 
smaller than the correlation length £(T), otherwise the 
magnetic configuration consists of many domains, and 
the bistability loses its meaning. We seek to identify the 
lifetime t of a magnetic state before it is flipped to the 
opposite magnetization direction. The flip is caused by 
thermal effects only and we are interested in a "loss of 
memory" , therefore we set B cxt = 0. It is expected, and 
found in the simulation, that r(T) follows to a good ap- 
proximation an Arrhenius law, 



r(T) =7t)fflq>(£k/fcBT), 



(4) 



where Eb is interpreted as an activation energy barrier 
and 1/to as an attempt frequency; fee is the Boltzmann 
constant (details on the method of calculation of r are 
given in Appendix A). Wc find that tq shows a devia- 
tion from a constant only at high temperatures, where 
the magnetization does not spend any appreciable time 
at a state before being reversed. What is important, 
however, is the dependence of Eb and To on the proper- 
ties of the chain. We keep a minimal model, considering 
only nearest-neighbour exchange coupling J = Jij and 
in most cases the same anisotropy K at each site (ex- 
cept when considering edge- anisotropy effects), whence 
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the energy- and time-scales of the system are defined by 
J, K, and A. Then, a domain wall that is formed in the 
chain will have a width of approximately £dw = J/K 
and a formation ener gy o f -Edw = 2\2JK (estimated at 
the continuum limit [29( of Eq. [T]). As we explain be- 
low, our main findings are categorized in the following 
regimes: (i) L <C £dw: the chain then behaves like a 
monodomain particle, as laid out already in the theory 
of Stoner, Wohlfarth and Brown [ll], [l2j], with the re- 
versal being controlled by a barrier of Eb = LK. (ii) 
L 3> Ldw: The magnetization reversal occurs due to a 
domain nucleation and subsequent domain wall propa- 
gation; the barrier entering Eq. (gj) is then Eb = Euw- 
Regime (ii) becomes relevant also for shorter chains as 
the anisotropy increases, especially above the so-called 
Ising limit K/J > 1, when the domain wall becomes 
atomically sharp; then even very short chains of a few 
atoms in length behave as in regime (ii). These results 
are presented in Section [TTJ Furthermore, in Section ITO1 
we investigate the possibility of different anisotropy pa- 
rameters in all three axes (i.e., we include also an in-plane 
anisotropy); in this case, entropic arguments imply that 
the attempt rate 1/Yo can be strongly affected, and this 
is verified by our simulations. In Section IIVI we discuss 
our findings for realistic system parameters of Co chains 
on Pt. We conclude with a summary in Section fVl 



II. CASE OF UNIAXIAL ANISOTROPY 

In the case of surface-supported chains the easy axis 
of magnetization, driven by spin-orbit coupling, is cither 
perpendicular to the chain or along the chain axis, as 
can be shown by symmetry arguments [30l . l3lj . For in- 
stance, in the case of Co chains on Pt [l| or Co islands 
and nanoclusters on Pt [32], [33| the easy axis is perpen- 
dicular to the surface (except in the case of chains at step 
edges, where there is a canting induced by the step edge 
environment ) 0, [H The anisotropy barr ler is 

typically between 1 and 10 meV, strongly dependent on 
the cluster size or shape. In this section we assume that 
the anisotropy tensor possesses rotational symmetry in 
the plane perpendicular to the easy axis; this particular 
Ansatz is fulfilled in free-standing chains but is not always 
valid when the chain is adsorbed on a surface. Thus, its 
approximate validity is strongly system-dependent and 
a more general case is investigated in Section [TTT1 As a 
rule of thumb, it should be more valid for chains that 
are contained in the first surface layer or for chains of 
biatomic cross section, compared to chains that are on 
top of the surface layer, because in the former case the 
in-plane asymmetry of the electron potential is smaller. 
Other effects that we neglect here are the possibility of 
temperature-dependent orientation of the easy axis that 
has been observed in certain systems H?} and a possi- 
ble canting of the anisotropy axis for the edge atoms, 
which has been found experimentally and theoretically 
[U [30L l3~il - l3"6j . However, at the end of the section we dis- 
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FIG. 1: (colour online) Simulation fits to the activation bar- 
rier height as a function of chain length taking into account 
a uniaxial anisotropy constant K. Cases K/J = 0.01 and 
K/J — 0.1 are shown. The horizontal lines correspond 
to the analytical solution of the continuum approximation 
E B /J = 2s/2^Kjj and agree rather well with the simula- 
tion results. For K/J — 0.01 we see a slight decrease of Eb at 
small lengths, indicating the onset of a transition to a Stoner- 
Wohlfarth behaviour as the length becomes comparable to the 
domain wall width. 

cuss some effects that could occur due to increased edge 
anisotropy. 

We proceed with a presentation of our results in the 
case of uniaxial anisotropy. We simulated the dynamics 
of monoatomic chains of various lengths up to a maxi- 
mum of 100 atoms, with K/J = 0.01 and 0.1; we chose a 
damping of A = 0.1. Values of K/J > 0.1 are rather ex- 
treme for supported transition-metal chains, but can be 
realized for example in cases of molecular magnets (where 
the exchange is comparably week), or for unsupported 
chains of 4d transition metal atoms [38j . The running 
time of the simulations varied, depending on the system, 
so that statistically reliable results could be obtained for 
the lifetime r. 

By fitting the Arrhenius law (jj) to our simulation re- 
sults we extracted Eb and verified the hypothesis that 
the energy barrier results basically from the formation of 
a domain wall, if the chain is not too short. The results 
presented in Fig. [1] show that Eb = £dw holds up to the 
limit L rj Low! for smaller lengths, a deviation is seen 
(in particular in the case K/J = 0.01, where Ldw = 20). 

Next we present in Fig. [5] simulation data on a chain 
of 100 spins for K/J = 0.1 at a temperature of fceT/J = 
0.11. The abscissa denotes the simulation time broken 
up in three time intervals; for each interval two panels 
are shown. In the upper panel, the ordinate shows the 
chain total magnetic moment along the easy axis (z-axis), 
M' ot , as a function of time. Reversal events are evident, 
indicated as B, D and E, when the total moment is re- 
versed. Each reversal event is relatively fast, followed by 
a longer time interval where the total moment is approx- 
imately stable at M* ot = ±85. However, we also see at 
A, C, and F attempted events that failed, with the to- 
tal moment returning to its previous state. In the lower 
panel of Fig. [5] we see a series of snapshots of the system 
over the same time span, visualizing the time-evolution 
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FIG. 2: (colour online) Time evolution of simulation snapshots for a 100-atom chain with anisotropy K/J = 0.1. The total 
simulation time interval is broken up in three sub-intervals of duration 0.25 x 10 6 h/J each. For each sub-interval, the upper 
panel shows the z-component of the total magnetic moment of the chain, Ml ot (t), while the lower, colour-shaded part shows 
a snapshot of the z-component of the magnetization at any site in time: green for Mi- Z (t) > 0, blue for Mi- Z {i) < 0. In the 
upper panels the attempts for magnetization reversal are witnessed by a significant change of (t), while in the lower panels 
one can see where the attempt is intitiated and how the magnetization-reversed region evolves in time. The letters and arrows 
correspond to events that are described in the main part of the text fSec. Hi)) . 



of the site-dependent magnetization. Here, the ordinate 
represents the position of the spin in the chain, i, while 
a colour-code is used for the sign of the local M i;z : green 
for Mi- Z > 0, blue for M^ z < 0. We can see how reversal 
attempts start by nucleation of a spin-fliped region that 
either propagates throughout the chain to end in rever- 
sal, as in B, D, and E, or is annihilated in an unsuccessful 
attempt as in A, C, and F. We also see that in some at- 
tempts the nucleation starts at the chain edge, as in B, C, 
E and F, while for others it starts in the chain interior, 
as in A and D. Although here we show a time- interval 
containing both kinds of nucleation, at not too high tem- 
peratures nucleation is more likely to occur at the edges 
than in the interior, because in the latter case two domain 
walls must appear simultaneously requiring the system to 
overcome a twice as high energy barrier (the latter pro- 
cess is known as soliton-antisoliton creation and has been 
studied e.g. by Braun [HI). For the same reason it is sta- 
tistically rare for more than one nucleation regions to be 
present simultaneously, at least as long as L -C £, there- 
fore we consider in our qualitative analysis only single 
nucleation events at the edges. Our arguments are sup- 
ported by calculations making use of a transfer-matrix 
approach to calculate thermodynamic properties of finite 
anisotropic Heisenberg [H showing that the edge-atoms 
of chains are subject to significantly stronger fluctuations 



than the atoms in the interior. 

The reversal time itself, i.e. the time that it takes for 
the domain wall to propagate through the system once it 
is created, is in general much shorter than the magnetic- 
state lifetime r. This effect is also seen experimentally 
by spin-polarized scanning tunneling microscopy of fer- 
romagnetic nano-islands [la]. There, fast reversal oc- 
curs during the scanning process, resulting in an image 
of seemingly two domains, while in reality the island is 
in a mono-domain state. 

The domain wall that occurs during nucleation per- 
forms a random walk back-and-forth in the chain, and 
eventually either returns to its starting point, resulting in 
an unsuccesful attempt, or reaches the other side, result- 
ing in reversal. The probability for the latter to happen 
in a random- walk model is [39[ p(L) = 1/L. Therefore, at 
constant temperature, the prefactor To of Eq. Q rises lin- 
early with L: tq(L) = aL, at least as long as Ldw *C L; 
this has been found also in simulations of Ising-modcl 
dynamics [40j |. More accurately, going beyond the Ising 
model, we must account for the fact that a domain wall 
centered at position x is attracted to the chain edge if it 
is too close, i.e., if x < |£dw or L — x < ^£uw- Thus the 
domain wall propagation corresponds to a random walk 
only between jLdw an d L — ^Ldw, so that we should 
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FIG. 3: (colour online) Magnetization lifetime r as a function 
of chain length L for the cases K/J = 0.01 and Kj J = 0.1 
(the latter was scaled up by a factor 10 for better legibility). 
The linear increase of r with L is due to the random- walk-type 
propagation of the domain wall after nucleation. The offset 
of the linear behaviour is due to the finite domain wall width. 
Vertical lines indicate the domain wall widths in the two cases 
(Ldw = 6 for K/J = 0.1, L DW = 20 for K/J = 0.01). 

observe an offset in the linear behaviour: 

t (L) oc (L - £ D w), for L > L DW . (5) 

This rule is reproduced by the results of our simulations, 
as shown in Fig. [3] where the calculated lifetime is plot- 
ted against the chain length. The offset in the 1st curve 
(K/J = 0.01) can be seen at lengths smaller or com- 
parable to the domain wall width (up to 20-25 atoms). 
Then, for L w LdWj the switching mechanism is interme- 
diate between nucleation and Stoner-Wohlfarth-like co- 
herent rotation, and a kink in the curve t(L) is seen. 
The offset is much reduced in the 2nd curve of Fig. [3] 
where K/J = 0.1, because the increased anisotropy leads 
to smaller Ldw ~ 6. A deviation from the rule 
in particular a saturation, is expected to begin as the 
chain length becomes too large approaching the correla- 
tion length. Then simultaneous occurence of more that 
one nucleation becomes statistically more and more rel- 
evant, until eventually the space-integrated magnetiza- 
tion vanishes at all times. Note that this linear increase 
([5]) of the lifetime with system size is in marked differ- 
ence to the behaviour of a Stoner-Wohlfarth-type three- 
dimensional particle, where the system size (number of 
atoms n) enters multiplicatively in the Arrhcnus barrier, 
Eb = nK, causing an exponential dependence of r on n. 
On the other hand, in the regime of validity of Eq. © 
(L > Ldw), the Arrhenius law (j4]) is still valid, but with 
Eb a constant (i.e., independent of the number of atoms). 

We close this section by an observation concerning the 
anisotropy of the edge atoms of the chain. It is well- 
known that, due to the different environment (only one 
neighbour) the anisotropy tensor of edge atoms, K cc i gc , 
can differ from the one in the chain interior, in both 
the direction of the principal axes and in the anisotropy 
strength, if e dge- These effects can affect the nucleation 
barrier as well as the nucleation frequency. However, 
there can be in principle an additional effect. Under "nor- 



mal" conditions a domain wall close to the edge is further 
attracted toward the edge and in the end expelled form 
the chain and annihilated, but if if edge becomes large 
enough, it will act as a barrier against the annihilation. 
Then the domain wall will be trapped in the interior of 
the chain for a longer time, eventually overcoming the 
barrier due to thermal fluctuations and becoming anni- 
hilated. Naturally, in such a situation the nucleation is 
also more difficult to achieve. We derived analytically an 
approximate condition for such an edge barrier by con- 
sidering the energy as a function of rotation of the edge 
spin and treating the rest of the chain in the continuum 
limit (see Appendix B). The condition reads: 

K cdgc > K + Ldw/4 (6) 

and is approximately reproduced by simulations that we 
performed. 

III. CASE OF TRIAXIAL ANISOTROPY 

As we commented earlier, the magnetic anisotropy of 
a chain on a surface is never really uniaxial, due to the 
asymmetry of the environment. We take this into account 
in the simplest way by introducing three eigenvalues of 
the anisotropy tensor, 

/ K x \ 
K = K y , (7) 
\ K z ) 

corresponding so to say to an easy, medium, and hard 
axis, so that the associated energy expression becomes 
— M • KM = -K x Ml - KyMl - K Z M%; the site index i 
has been suppressed here, and it is assumed that all sites 
have the same anisotropy. Without loss of generality we 
fix K z > K y > K x = (the latter equality can be ad- 
justed by a shift of the energy zero that does not affect 
the equations of motion, while the order of the inequali- 
ties will not affect the statistical results on the lifetime, 
merely the path that the spins follow to flip). Thus, in 
effect the two energy minima of the single-ion anisotropy 
are along the ±z-axis, separated by a low anisotropy bar- 
rier Ai ow = (K z — K y ) if the spin rotates in the y-z-plane 
or by a high anisotropy barrier Ahigh = K z if the spin 
rotates in the x-z-plane; effects of the exchange J come 
on top of this as before. 

The question that we seek to answer is how the pres- 
ence of the second anisotropy barrier affects the mag- 
netization lifetime. To make the question more concise, 
we wish to compare the lifetime in the previously cal- 
culated uniaxial case, where Ahigh = Ai ow = K, to the 
case Ahigh > Ai ow = K. The minimal barrier for domain 
wall formation in long chains or for coherent rotation 
in short chains remains 2\/2JK or LK respectively, so 
that the exponent Eb in the Arrhenius law (j4]) should 
not be affected. However, the attempt frequency t^ 1 is 
expected to change, as the energy lanscape of the paths 
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that connect the two minima is now different. For exam- 
ple, in the case of uniaxial anisotropy both Bloch-type 
and Neel-type walls are equal in energy, while in the tri- 
axial case one of the two types will correspond to a higher 
energy; this must give a difference in entropy. 

As it turns out, for a given Ai ow , the lifetime decreases 
with increasing Ahigh; our findings are summarized in 
Fig. |U This is at a first glance counter-intuitive: the 
fluctuating force enforces a random walk, regulated by 
the barriers, that has to climb up from an energy mini- 
mum to a maximum or saddle point before descending to 
the opposite minimum. Increasing Ahigh means that part 
of the lanscape is visited with smaller frequency, so that 
there are less escape paths. Then part of the random- 
walk steps would be, so to say, waisted in attempts to 
climb up the high barrier. However, contrary to the case 
of a confined particle, the latter argument is not valid in 
the case of spins, because of the precessional motion. The 
spin trajectory has a component following the direction 
of the effective field (toward the energy minimum) , given 
by the damping term of Eq. ([2]), but the stronger, pre- 
cessional part drives the spin perpendicular to the energy 
gradient. Therefore, whenever the fluctuating field forces 
the spin in the direction of the high barrier, the precession 
drives it in the direction of the low barrier. In this way, 
the random-walk steps in the direction of the high barrier 
are not waisted but converted into random- walk steps in 
the direction of the low barrier. The shape of the result- 
ing trajectory approximates an ellipsis with its long axis 
along the low-barrier direction. The random walk be- 
comes in a sense more and more one-dimensional, as the 
difference Ahigh — Ai ow increases, and it is known that 
random walks in one dimension propagate faster away 
from the starting point than in two dimensions. There- 
fore the magnetization lifetime becomes shorter. In the 
limit Ahi g h 3> Ai ow , the situation resembles that of an 
y-z-easy-plane model with a hard axis in the x direction, 
which obeys the dynamics of precessional switching, long 
known in micromagnetics (4lj . At this limit, even a small 
deviation of the spin toward the x direction due to the 
fluctuating force will cause a lengthy precession around 
x, in the y-z plane, which can even last a few rotations 
if the damping is very low. 

The behaviour is visualized for a single spin in Fig. [5] 
On the left panel we see a top-view of the trajectory of 
a spin under a uniaxial anisotropy, with K x = K y = 0, 
K z = K > 0, starting from the state M x = M y = 0, 
M z = 1. The spin is driven away from equilibrium by 
a fluctuating field, and the trajectory is, on the aver- 
age, isotropic in x and y. On the right panel we see the 
trajectory for the same simulation time and under the 
same conditions, with the exception that here A[ ow = K , 
Ahigh = 10K. The point is not just that the trajectory 
is elliptical instead of isotropic, but also that in the same 
time interval the spin has managed to reach closer to 
a "turning point" of the barrier, where + My = 1, 
M z = 0. 

In the case of chains (Fig. Ql we observe in the top 
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FIG. 4: (colour online) Magnetization lifetime as a function 
of temperature for a 100-atom chain for several single-ion 
anisotropy barriers, including triaxial anisotropy (top, linear 
scale; bottom: logarithm of the lifetime vs. inverse tempera- 
ture). The low barrier, Ai ow = 0.01J, is kept constant, while 
Ahigh is varied between Ai ow and 10Ai ow . A significant reduc- 
tion of the lifetime is observed even at a moderate difference 
of 10% between Ahigh and Ai ow . At very high difference a 
saturation is seen. The activation barrier is in all cases the 
same, equal to the domain wall energy, Eb = 2\/2JAi ow , as 
can be seen in the lower panel where all data sets show the 
same slope. 




FIG. 5: (colour online) Trajectory of a single spin under 
anisotropy and fluctuating force. Left: Uniaxial anisotropy, 
K x = K y = 0, K z = K > (Ahi g h = Ai ow = K). Right: 
K x = 0, K y = 9K, K z = 10JC (Ahi g h = IQK, A low = K). 
Both simulations correspond to the same total simulation 
time starting from the state M x — M y = 0. Evidently, in 
the presence of different barriers in x and y directions (right 
panel) the spin reaches faster the turning point at the edge of 
1, M z = 0). 
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panel a drop of the lifetime with increasing Ahi g h, while 
keeping Ai ow = 0.01J fixed. The curves show the life- 
time as a function of temperature. At all temperatures, 
the longest lifetime is observed for Ahigh = Ai ow , and it 
is evident that even a small increase of 10% in Ahigh can 
affect the lifetime quite noticeably. When Ahigh is twice 
the value of Ai ow , the lifetime is reduced by an order of 
magnitude. However there appears also to be a satura- 
tion: further increase of Ahigh to nine times the value 
of Ai ow makes no significant difference any more for the 
lifetime. 

In the lower panel, where the same simulation data are 
shown in a log tvs. 1/T plot, the Arrhenius behaviour as 
well as the dependence on Ahigh — Ai ow becomes easier to 
see also at higher temperatures; evidently the activation 
barrier is the same in all cases (and equal to 2-^/2 JAi ow ), 
as can be seen from the fact that the slope of log tvs. 1/T 
does not change. Note that at low temperatures the data 
have larger error bars due to the fewer switches that occur 
in a given simulation time interval. 



IV. REALISTIC PARAMETERS AND RESULTS 
FOR Co CHAINS ON Pt 

Now we relate our general analysis to realistic sys- 
tems. For the typical system of supported monoatomic 
Co chains and nanoislands, for example, values of J of the 
order of 20-60 meV on Cu, Au and Pt surfaces have been 
calculated and reported in the literature 0, l42Tl44l | , de- 
pending on structural parameters such as island or chain 
geometry, size and substrate. Similar is the case with 
Fe nanostructures. The anisotropy, on the other hand, 
is found in experiments to vary by two orders of magni- 
tude depending on structural parameters [3^ , from the 
sub-meV range for large clusters up to values as high as 
10 meV for single adatoms; for monoatomic Co chains on 
Pt(997) step edges, a value of 2 meV has been found by 
experiments [l|, while the value drops to 0.34 meV for 
biatomic chains. Clearly, the anisotropy is much more 
sensitive than the exchange coupling. Concerning the 
appearance of two anisotropy barriers, Ai ow and Ahigh, 
calculations [Ifl on Co monoatomic chains on Pt(lll) 
yield Ai ow = 0.8 meV and Ahigh = 1 meV; this 25% 
difference between Ai ow and Ahigh can play a significant 
role on the lifetime, as our calculations of Sec. Illll show. 

Following the above observations we are in a position 
to make an order-of-magnitude prediction for the life- 
time of Co monoatomic chains on Pt. We assume val- 
ues of J = 50 meV, A low /J = 0.016 A hig h/J = 0.02. 
The only parameter that is at this point arbitrary is the 
damping, here chosen to be A = 0.1, which is not ac- 
tually known for the case of Co monoatomic chains on 
Pt. In permalloy, for example, the damping is one or- 
der of magnitude smaller J46(. However, Pt is charac- 
terized by a much stronger spin-orbit coupling, which 
can enhance the damping (for example, FePt nanoparti- 
cles are reported to have a damping parameter as high 



as 0.76 Hll)- Our current choice has to be corrected 
once there is more experimental or theoretical data, in- 
cluding exchange parameters beyond first neighbours and 
an account of the anisotropic exchange (Dzyaloshinskii- 
Moriya interaction), on the particular type of systems. 
Keeping this in mind, we find for a 50-atom chain an 
inverse attempt rate of tq ~ 0.6 ps and a barrier of 
Eb = 2y / 2JA\ ow pa 18 meV to be substituted in the 
Arrhenius law ((4]) . Assuming detection of the reversal by 
a scanning tunneling microscope (see e.g. Refs. [l5| and 
j48| ) , and given that this method has a relatively low time 
resolution of the order of msec [49| , we consider a lifetime 
of 1 second as reasonable for a reliable experimental anal- 
ysis of spontaneous reversal. Then our simulation sug- 
gests that experiments on 50-atom-long chains require 
temperatures of approximately 10 Kelvin or lower. 



V. CONCLUSIONS 

We have performed classical simulations of atom- 
istic spin dynamics for supported monoatomic magnetic 
chains, with focus on the magnetization lifetime between 
temperature- induced subsequent magnetization reversals 
at chain lengths below the correlation length. The cal- 
culations were carried out by integrating the Landau- 
Lifshitz equation of motion of classical spin systems in 
the presence of stochastic forces. 

In summary, the main findings of this work are the fol- 
lowing, (i) For long chains (longer than the width of a 
domain wall), the lifetime is governed by an Arrhenius 
law @, with a barrier equal to the domain wall energy. 
For short chains, the barrier eventually becomes equal to 
the anisotropy energy, (ii) For long chains, the mecha- 
nism for reversal is governed by nucleation of a region of 
reversed magnetization at the chain boundary, forming a 
domain wall, and propagation of the domain wall through 
the chain via a random-walk-like procedure. This yields 
a linear dependence of the lifetime on the chain length 
for chains longer than the domain-wall width, contrary 
to an exponential dependence that would be expected for 
a Stoner-Wohlfarth type of system, (iii) In case of a tri- 
axial anisotropy tensor with non-degenerate eigenvalues, 
i.e., with a low and a high anisotropy barrier for magne- 
tization reversal, the lifetime is reduced compared to the 
case of only a low barrier. This effect is related to the 
precessional motion of the spin, which tends to transform 
fluctuations in the high-barrier direction into fluctuations 
in the low-barrier direction. A triaxial anisotropy should 
always be present in surface-supported chains, as can be 
argued on the grounds of symmetry of the system. 

As far as fundamental science is concerned, the most 
important question to be answered by this type of calcu- 
lations is perhaps the one of the limitations in size, time 
and temperature, that the classical spin dynamics is still 
a valid approximation, especially in metallic systems. In 
principle, all ingredients that enter the Landau-Lifshitz 
equations are experimentally accessible; even more so for 
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chains, where the geometry is simple enough so as not to 
introduce further arbitrariness. The result of the simu- 
lations, in particular the magnetization lifetime, is also 
accessible by experiment. A systematic study in this di- 
rection will be highly interesting. 
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Appendix A: Some details on the method of 
calculation 

Equation @| is a stochastic differential equation 
with multiplicative noise which is interpreted in the 
Stratonovich sense. The fluctuating force fi(t) has 
white-noise properties, i.e. (i) its time-average van- 
ishes: (fi) = 0; (ii) it is fluctuating much faster 
than the characteristic precession time so that time- 
separated and space-separated fluctuations are decorre- 
latcd: {[fi(t)] a [fj{t')b) = 5ij5 a(3 6(t ~ t'), where a,/3 = 
x,y,z denote the vector components and i,j are atomic 
sites; and (iii) its amplitude e is connected to the damp- 
ing via the fluctuation-dissipation theorem: e 2 = 2Xk^T 
with &b the Boltzmann constant. Eq. ^ is integrated 
via a weak Runge-Kutta method which was suggested by 
Milstein and Tretyakov [50| . The error scales with the 
step size as Oih^ + e 2 h 2 ). Further information can be 
found in Ref. [24j . 

In practice, at each Runge-Kutta step and for each spin 
at site i, a fluctuating field fi is calculated by a set of 
random numbers obeying an appropriate distribution so 
that the above requirements (i— iii) are fulfilled. Require- 
ment (iii), in particular, guarantees that the fluctuating 
field intensity corresponds to the particular simulation 
temperature T. The field fi acts in addition to the ef- 
fective field Bf and to the damping term as shown in 
Eq. (|2|). The random nature of fi comes into notice in 
the noise observed in the time-dependent magnetiation in 
Fig. [5] For the calculation of r(T), during a sufficiently 
long simulation time at a fixed temperature, successive 
time intervals Ati between N successive reversal events 



(as the ones seen in Fig. [2]) are recorded, and at the end 
averaged as t(T) = $3 i=1 Atj/iV. 

Appendix B: Justification of Expression ([6]) 

Consider a domain wall where the magnetic moment 
of each atom i has an angle 8i with respect to the easy 
axis (the easy axis is taken to be the same for all atoms). 
Suppose that the wall is not at the interior of the chain, 
but close to the edge, and that it is formed by constrain- 
ing the moment of the first atom (labelled "0") by 9o; 
we accept the boundary condition at the deep interior of 
the chain that 0i = for i — > oo. Then, a boundary con- 
dition at the edge 9q = n corresponds to a fully formed 
domain wall, while 9q = corresponds to the ferromag- 
netic ground state. 

Now we take a fixed #o, < 9q < ir. The result- 
ing structure will have the form of a domain wall that 
has been abruptly cut at i = 0. If we assume in the 
discrete model a nearest neighbour distance of length a 
and employ the continuum approximation p9j . then the 
site index i changes to a continuous variable x/a and we 
obtain an explicit energy functional for a chain starting 
with a restricted angle of #o at x = 0: 

E(9 o) = f d x (±Ja (£) 2 + f sin^))). 

(Bl) 

Minimizing the energy functional under the boundary 
conditions 9(0) = 9 and 9(oo) = leads to an energy of 
the domain wall of 

£dw(0o) = V2JK(1 - cos0 o ), (B2) 

which will be used as an approximation to the discrete 
model. 

Next we increase the anisotropy of only the edge atom 
from K to -Ked ge - This change affects only the energy of 
the edge atom, changing it by — i4f e dg e cos 2 9 + K cos 2 9 , 
while the configuration remains the same. We have now 
an energy 

E(6 Q ) = V2JK(1 - cos 6> ) - K edge cos 2 6 + K cos 2 9 . 

(B3) 

For the domain wall to be trapped in the interior of the 
chain, this expression must have a maximum for < 9 < 
ir, corresponding to the high-point of an energy barrier 
upon rotation of the edge spin. Differentiation of (|B3[) 
with respect to $o shows that such a maximum exists 
under the condition 

K > y/JK/2, (B4) 
which is equivalent to Expression (|6]). 
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